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Abstract 

We investigate the influence of interference between Coulomb interaction and 
impurity scattering on the static electronic response x(0, q) in disordered met- 
als to leading order in the effective Coulomb interaction. When the transport 
relaxation time r^^ is much shorter than the quasiparticle life time, we flnd a 
sgn(2pp — q)l \/\2pp — q\ divergence of the polarization function at the Fermi 
surface {q = 2pp). It causes a phase shift of the Friedel oscillations as well as 
an enhancement of their amplitude. Our results are consistent with experi- 
ments and may be relevant for understanding the stability of the amorphous 

state of certain alloys against crystallization. 
PACS numbers: 61.25.Mv, 64. 70. Kb, 71.45. Gm, 71.45. LR, 72.15.Cz 
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Due to the restriction of electrons to states with momentum less than or equal to the 
Fermi momentum pp^ electronic density correlations exhibit spatial oscillations, generally 
referred to as Friedel oscillations [1]. It has been suggested [2-4] that the resulting oscilla- 
tory contribution $°*^(r) to the effective pair potential between the ions in a metal could 
be important for stabilizing the amorphous state of certain alloys, in that the ions arrange 
themselves in commensurate positions at the minima of $°*^(r). Indeed, experiments on a 
wide class of quenched noble-polyvalent metal alloys, which exhibit a crystalline to amor- 
phous transformation as a function of noble metal concentration, clearly show a relation 
between atomic and electronic structure [4,5]: In all those systems the transition occurs, 
when the Fermi energy ep lies in a structure-induced pseudogap of the density of states 
(DOS). Moreover, in the transition region, the atomic spacing a coincides with the Friedel 
wave length Xp = tt/pf over distances of several Xp. However, the atomic positions are 
shifted by Xp/A compared to the positions of the minima of $°*^(r) as calculated for a clean 
electron sea. Assuming that in the amorphous state the atoms are sitting in the minima 
of $°*^(r), it may be conjectured that in those systems Friedel oscillations acquire a phase 
shift. The problem of the strength of the Friedel coupling has been a subject of vivid in- 
terest. In disordered systems the oscillations of the impurity-averaged potential ($°*^(r)) 
are exponentially damped [6] due to a random phase, which arises from the spread of the 
Fermi wave length over a width given by the elastic mean free path [7,8]. Yet, Jagannathan 
et.al. [9] showed that higher moments of the distribution of $°*^(r) produce a long-range 
coupling. In this way, the behavior of RKKY spin glass systems [10] at low temperatures T 
could be explained. In the amorphous systems considered here, however, thermal energies 
(r = 77 K) are typically of the same size as the depth of the first minimum of $°*^(r) of 
even a clean electron sea. Thus, the binding of ions by $°*^(r) and the existence of a phase 
shift of Friedel oscillations have remained controversial issues. 

In this paper we show that the interplay between Coulomb interaction and impurity 
scattering can give rise to a phase shift as well as an enhancement of the amplitude of 
the electron density oscillations in systems where the electronic transport time Tfr is short 



compared to the quasiparticle life time r^p. We give arguments that this condition is sat- 
isfied near the crystaUine to amorphous transition. The above efi'ect is intimately related 
to a previously undiscovered instability of a disordered electron sea towards formation of 
concentrical charge density oscillations in 3 dimensions. 

While the response of disordered electron systems has been studied extensively in the 
hydrodynamic limit [11], Friedel oscillations are governed by the static electronic response 
x(0,(j') at momenta q ~ 2pi?, which is unexplored to a large extent. In a disordered Fermi 
liquid the interaction vertex is dressed by impurity ladder vertex corrections (Fig. 1(a)) 
[12]. Due to their diffusion pole structure, any two-particle interaction acquires a strongly 
retarded, long-range part [13,14,11]. This can drastically alter the screening as shown 
below. If one neglects any coherent scattering processes, impurity vertex corrections are 
simply inserted into the vertex of the well-known random phase approximation (RPA). 
When taken into account on this classical level of approximation, neither substitutional nor 
structural (amorphous) disorder produces a phase shift of the density oscillations [15]. The 
reason is that in those diagrams the impurity ladder vertex F enters in the static limit 
at finite q, so that the diffusion pole is suppressed. We are thus led to consider vertex 
corrections where F is integrated over frequency and momentum, i.e. does give an infrared 
divergent contribution. There the hydrodynamic transport properties become important, 
although the response is taken at large external momenta. 

The strong isotropic density correlations between the ions of the system [5] have impor- 
tant implications on the electronic transport: The impurity scattering T matrix tppi has a 
pronounced peak at a momentum transfer of size q = \p' — p\ = qo, where qo := 2tt j a. This 
can be seen directly from the ionic structure factor S{q), measured by electron diffraction 
[5], which is proportional to tpp'. The scattering angle for "on-shell" scattering (i.e. in 
and outgoing momenta p, p' on the Fermi surface) with momentum transfer qo is given by 
COS0 = 1 — 2{qol2ppY . It is seen that backscattering is enhanced when qo ~ 2pp, or a = A^?, 
in analogy to Bragg refiection in 1 dimension. The transport rate r^^ is obtained by solving 



[16] the impurity ladder Bethe-Salpeter equation as Tj-^^ /t J- = 1 + TqpM, where M denotes 
the "current kernel", 

Gp{0~) and Vp are the advanced Green function at the Fermi energy and the group velocity, 
respectively. A " denotes the imaginary part throughout this paper. We assume |tp^p'P oc 
^{\p' ~ P\ ~ <?o)- Because of the peak structure of G'' at the Fermi surface, M is governed 
by on-shell scattering, and in the transition region to the amorphous state (qo ~ 2pi?) r^^^ 
is enhanced compared to t~^ because of the antisymmetry of the current vertices (vpq). To 
obtain an estimate of Tj7 /t~^ we have evaluated Fq. (1), where the single-particle quantities 
have been calculated in self-consistent Born approximation. Note that this constitutes a 
conserving approximation. It reproduces the observed [5] structure-induced pseudogap of 
the DOS (Fig. 2). Fig. 2 also shows a strong enhancement of t^^^ /t~J- when 2pF — qo, or 
equivalently, when ep is near the minimum of the DOS. This will become important below. 
We now turn to calculating the polarization function I[[(jj^q). All corrections to I[[(jj^q) 
to leading order in the effective Coulomb interaction Vq[z^Z) are shown in Fig. 1(b). We 
note that these diagrams are generated in a conserving approximation from the same class of 
free energy diagrams that produces, on the single-particle level, the Al'tshuler-Aronov [13] 



self-energy, which causes a y\E — cp] singularity of the DOS in a disordered Fermi liquid, 
as was first shown by those authors. n(c<;, q) is related to the density response x(^7 <?) t»y 
X = n/(l — Vqll)^ with Ug, the bare Coulomb interaction. For complex frequencies z -\- Z 
and z of in and outgoing lines, Vq{z^ Z) is given by 

V Attc 

where e^^'^^Z^q) = 1 + 27rz a sgnZ"/{Z + iq^D sgnZ") is the disordered RPA dynamical 
dielectric function and 



T{z,Z,q) 



iJTtr sgnz" 

„2 7 



if z"{z + Z)" < 



Z+ig^D SgnZ" ^ ' /o\ 

1 otherwise 



is the dressed interaction vertex, D = 1/3 v^Ttr and a = ne'^Ttr/m denoting the diffusion 
coefficient and the Drude conductivity, respectiveiy. In Eq. (2) e^^^ and F^ combine to a 
single diffusion pole indicating the retarded, long-range character of Vq. Since the Drude 
conductivity, expressed in terms of the plasma frequency cOp as a = cuVtr, is typically 
three orders of magnitude larger than the Fermi energy e^?, the nondivergent terms may be 
neglected. 

It may be shown by Fourier transformation that a phase shift of the Friedel oscillations 
corresponds to a peak structure of xlO?*?) 3.t q = 2pp. As seen below, the appearance of 
such a structure is connected with a divergence of the effective interaction at zero frequency 
and momentum transfer from the particle to the hole line in 11(0, (j'). Since this divergence 
is obviously strongest in diagram (1), Fig. 1(b), we only need to consider this diagram, 11^^^ 
in order to calculate the leading singular behavior at the Fermi surface [q = 2pp). Note in 
this context that, while diagrams (l)-(3). Fig. 1(b) cancel each other exactly in the limit 
cu, (j' — )• [14], this is no longer the case for finite momenta. 11^^^ reads 

n(^)(0,?) = -2r2^^G'p+_(z^) Gp-4iE) Vq,{iE,iu') Gp++{iE + iu) Gp-+{iE + iu) (4) 

E,cu' p,g' 

^,J TT (ep+_ - ep__)(cu' + ep+_ - ep++)(cu' + ep+_ - ep_+) 
+ nondivergent terms. 

Here, / is the Fermi function, p^ stands for p-\-q/2 — q' /2 etc., and V indicates the principal 

value. The io' integral over the singular part of v in Eq. (4) extends up to frequencies 
~ (27r/3)Tj7^. Near the amorphous transition t~J- is substantially smaller than r^^^ (see 
above). Thus, in evaluating 11^^^ we have assumed the quasiparticle pole of Gp{z) to lie 
on the real axis for the sake of clarity of the present calculation. Careful inspection of all 
terms contributing to 11^^^ shows that only the one written explicitly in Eq. (5) is divergent 
at the Fermi surface. Due to the strong q' = divergence of v'', we can set q' = in all 
nonsingular terms [13]. The momentum integrals can then be performed without further 
approximations: Applying an integration by parts in the integral over the component of 



p parallel to q elucidates the fact that Il^^> is governed by contributions from the Fermi 
surface, and we obtain (T <^ ep) 



na)(o,,)= ^^ ' r^ '^ "^'^^ 



16 {2eFTt,yl^ J-e^ a:3/2cosh2^ 



2T 



Fi[u^ x) — F2{u^ x) 



I, ,15/2 

Fi(M,a;) = iv(a;) — (1 )L{xi), F2{u,x) = \/l — Axu \L{x) — (1 )L{xi)\ 



where iv(x) = n(°)(0,?) = -2mpF/i2TrY [1 + (1 -x'^)/{2x) ln|(l + x)/(l - x)\] is the Lind- 
hard function {T,uj = 0), and x = x{u) = {q/2pF)/ ^Jl + u/ep, Xi = x — u, x = x/y^l — Axu 



and Xi = (x — u)/\/l — Axu. Each term in the first integral over u has a nonintegrable diver- 
gence at M = 0, which is compensated only by a cancellation up to 0{u) when Fi and F2 are 
combined. It follows that the q dependence of Il^^> is governed by the second derivative of 
L{x), which diverges as (1 — x)~^ . The second u integral in Eq. (6) is nondivergent and may 
be dropped. To make the q dependence explicit, it is convenient to substitute u = \x — l\y. 
In this way II^^^ may be written in the vicinity of x = 1 as 

* '" 16(2tFT„)'/2 y cosh' Jy ^TTTj ' '' 

where 

.,,.„,„ ,, , \f.(ii..x\ - FJii..x\\ 

u=\x-l\y ^o\ 



ilix\x-i\ ^y [Fi{u,x) - F2{u,x) 



3 \y\^^'^ x-^/^ [x — 1) 

is a continuous, positive function at x = 1. The integral may be evaluated numerically as 
/(I) ~ 11.8{2mpF/{27rY). 

Eq. (7) is the main result of this paper. It is seen that at T = the static polarization 
function 11(0, q) = n^°^(0, q)-\-Il^^'{0, q) acquires a powerlaw divergence with a sign change at 
q = 2pp, the half integer exponent being characteristic for diffusive behavior. Consequently, 
below a temperature r< the static dielectric function e{q) = 1 — ^^11(0, q) becomes less than 
1 near q = 2pp, and it follows from a Kramers-Kronig relation that the dynamical stability 
criterion, c<;x"(c<; — zO, q) > 0, is no longer fulfilled for all frequencies: The disordered electron 
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sea becomes unstable towards formation of charge density oscillations with wave length 
Xp = t^/pf- Note however, that a prerequisite for this instability is r^^^ ^ t~^, which means 
in the present case that there are ionic correlations with wave number qo = 2pp. Hence, the 
negative 2pp peak in x(0, q) is compensated by a peak in the static ionic response, which is 
equal and opposite in sign due to local charge neutrality. Further calculations will be needed 
to confirm that the ionic and the electronic structures are indeed mutually stabilized in this 
way. 

Our results are summarized in Figs. 3. and 4. The screening charge density p{r) is 
computed as the Fourier transform of Ij e{q) — 1 [17]. In p{r) the disorder induced oscillations 
are superimposed on the conventional cos[2ppr — (^), (^ = 0, Friedel oscillations. As long as 
e(<?) ^ 1, the 2pp powerlaw singularity results in a steplike behavior of l/e(<?) at finite T, 
however with an opposite sign compared to the clean response function (Fig. 3(b)), and hence 
produces density oscillations oc — cos[2ppr). Thus, there is no significant phase shift (p in 
this region. Rather, the Friedel oscillations are weakened and even become overcompensated 
by the quantum corrections. In Fig. 4 the latter is seen as the point where the amplitude A 
vanishes and the phase shift jumps to (p = tt. As epTfr is further decreased, l/e(<?) develops 
an increasingly asymmetric peak structure at 2pp (Fig. 3(c)), which leads to a contribution 
oc sin[2ppr) superimposed on the conventional oscillations. Thus, as the weight of the 2pp 
peak grows, the amplitude A of the density oscillations increases and (p approaches 7r/2. 
Fventually, at a critical transport rate Tj^\, e{2pp) becomes 0, and the peak turns into a 
divergence, i.e. the amplitude of the disorder induced density oscillations diverges, marking 
the breakdown of perturbation theory. Near that point, the conventional Friedel oscillations 
are negligible, resulting in a phase shift oi p = 7r/2. It should be emphasized that, although 
the divergence of 11(0, (j'), Fq. (7), is changed to a peak when a finite quasiparticle life 
time Tqp is included, the divergence of l/e(<?) persists as long as r^^^ ^ t~^. From the 
above discussion it is now clear that p ~ 7r/2 is intimately related to a strong enhancement 
of the amplitude of the oscillations. This explains the curious experimental finding [4,5] 
mentioned in the introduction that the transition to the amorphous state occurs at a shift 
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of the atomic positions by Ai?/4, corresponding to (p = 7r/2, because at that same point 
the effective pair potential $°*^(r) becomes particularly strong and dominates the crystal 
potential and thermal energies. 

In conclusion, we have calculated the static electronic response of a disordered metal, 
including quantum interference between Coulomb interaction and impurity scattering to 
leading order in the effective Coulomb interaction. In a situation where the transport relax- 
ation rate is much greater than the quasiparticle relaxation rate, the quantum corrections 
display a previously unknown powerlaw divergence at q = 2pp. This structure leads to a 
phase shift as well as to an enhancement of the amplitude of the Friedel oscillations. Hence, 
this effect seems to explain experiments [4,5] on certain strongly disordered alloys, and, in 
particular, may act stabilizing on the amorphous state of these materials. A detailed ex- 
position of the calculations presented here, including finite quasiparticle life time will be 
presented elsewhere. 

We are grateful to P. Wolfie, A. A. Aronov, N. W. Ashcroft, and P. Haussler for valuable 
discussions and comments. This work was supported in part by SFB 195 (T.K.) and by the 
Alexander von Humboldt Foundation (J.K.). 
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FIGURES 
FIG. 1. (a) Impurity dressing of the interaction vertex, (b) All corrections to the polarization 

function Il(uj,q) in leading order in the effective Coulomb interaction Vg are shown. Vg (Eq. (2)) is 

denoted by a bold wavy line with solid triangles at the ends. The shaded blocks represent diffusion 

vertices. In diagrams (2)-(5) a corresponding contribution is obtained by interchanging the particle 

and the hole line. 

FIG. 2. T~J- , T^^ in units of g^/2m, and the DOS N(E) (arb. units) as a function of energy. A 
quadratic band E = p"^ /2m is assumed for the clean electron sea, and parameters are chosen such 

that T~^ ~ 0.12ep at the Fermi surface. : clean DOS N'-°''(E). The vertical line indicates 

the position of the Fermi level for 2pp = Qg. 

FIG. 3. l/e(q) is shown for a clean (a) and a strongly disordered (epTir = 1.8 (b) and 1.23 (c)) 
electron liquid as calculated from Eq. (7) near q/2pp = 1. All numerical evaluations are done with 
an electron gas density parameter [17] of r^ = 4 and at T/ep = 0.001. 

FIG. 4. Phase shift cp and amplitude A of the 1st maximum of p(r) as a function of ep(Tir — Tir^c)', 
epTir^c = 1.217. The inset shows p(r) (arb. units) as a function of distance r from an impurity ion 
(a) in the clean and (b) in the strongly disordered case (epTir = 1.23). The phase shift and the 
enhancement of the amplitude are clearly seen. 
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